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The intrinsic contribution to the spin Hall effect in a 2DEG with non-magnetic impurities is studied 
in a quantum Boltzmann approach. It is shown that if the steady state response is perturbative in 
the spin-orbit coupling parameter A, then the precession term — vital for Dyakonov-Perel relaxation 
and the key to the spin Hall effect in previous similar Boltzmann studies — must be left out to first 
order in spin-orbit coupling. In such a case one would have that to lowest order in the parameters 
electric field, spin-orbit coupling, impurity strength and impurity concentration there is no intrinsic 
contribution to the spin Hall effect, not only for a Rashba coupling but for a general spin-orbit 
coupling. To cover all possible lowest order terms we consider also electric field induced corrections 
to the collision integral in the Keldysh formalism. However, these corrections turn out to be of 
second order in A. For comparison we derive some familiar results in the case when the response is 
not assumed to be perturbative in A. We also include a detailed discussion of why a relaxation time 
approximation of the collision integral fails. Finally we make a comment on pseudospin currents in 
bilayer graphene. 



I. INTRODUCTION 

The idea in spintronics is to manipulate the electron spin for information storage and transfer, as done with the 
electron charge in electronics. Possible advantages could for example be smaller resistive losses, or the additional 
richness that comes from the non-scalar nature of spin. In the study of spin currents particular attention has been 
given to the spin Hall effect (SHE) in a two-dimensional degenerate electron gas (2DEG) with spin-orbit coupling. 
Perpendicular to an applied electric field, opposite spins travel in opposite directions, thus creating a spin current 
without a net charge current. Such a spin-current can lead to an accumulation of spin at the edges of a sample, 
although in contrast to the electrical charge analogy this is not given. An individual spin can change its direction, 
for example due to spin precession, and therefore the accumulated spin polarization is not a conserved quantity. The 
SHE opens up one possibility of manipulating spin with electric fields. 

The experiments on the SHE are few and recent. [l|, 0, H, H, H, Q In contrast, much work have been devoted to 
the theoretical aspects (see e.g. the reviews 0, H, @) and even in recent years there has been an intense discussion 
about the different mechanisms and how to compare results reached by different formalisms. This discussion is closely 
related to the one on the anomalous Hall effect ( AHE) [l3| , for example in the distinction between extrinsic (impurity 
related, e.g. skew scattering, side-jump etc) and intrinsic (band structure related) mechanisms. The latter could be 
due to structure inversion asymmetry of the confining potential (leading e.g. to a Rashba coupling) or due to bulk 
inversion asymmetry (allowing for a Dresselhaus coupling). 

The earliest theoretical works on the SHE were done by Dyakonov and Perel in 1971 [ll|, who showed that a 
spin-orbit coupling leads to a spin current perpendicular to the electric field. The name Spin Hall effect was coined 
in 1999.(13 A universal spin Hall conductivity in a 2DEG was proposed. [H, [ll, [l5(] However, impuritites were 
neglected in these studies. Many studies have later shown that universality is lost when impurities are taken into 
account. [11, [13, EES ill, 112, 111, 113, 111, [2^^ On the other hand, there can still be topological edge modes, 
responsible for the quantum SHE (see e.g. ref. 1291). 

The effect of non-magnetic impuritites on the intrinsic contribution to the SHE has been studied using a Boltzmann 
approach [H, [13, [H, [ll, [H, [SI as well as diagrammatic methods [H, [H, [H, [H, [H, [13, [H like the Kubo formalism. 
The diagrammatic methods are more systematic and have a more general range. The Boltzmann approach offers more 
intuition, but has typically been implemented by identifying distinct processes or hand-picking different contributions 
rather than systematically covering all possible contributions through a formal apporach. Within the Boltzmann 
language one finds several different approaches. We are going to use the quantum Boltzmann approach, which treats 
spin coherently. 
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This paper is going to deal with the intrinsic spin Hall effect in a steady-state calculation to first order in spin-orbit 
interaction. We attempt in a Keldysh derivation of the Boltzmann equation to account for all contributions present 
to lowest order in electric field, in spin-orbit splitting and in strength and concentration of non-magnetic impurities. 
(Sees. IllUfVl and sec. [X]) The spin-orbit interaction Hso = cr ■ b is chosen to be of the isotropic form b = b{k)b{9), 
where the unit vector b has a winding number A^, i.e. + iby = e^^a+tNe ^^^[^1^ being a constant. 

Many recent theoretical studies come to the conclusion that the spin Hall effect vanishes for N = ±1 (e.g. for 
the Rashba and linear Dresselhaus spin-orbit couplings) for point-like impurities [l^, [iM 0i HO, HH, [13, [2^, 
[2^, [23, [13 as well as finite range impurities [13, [3, [3, [ll| ■ General arg uments have been proposed to explain this 
vanishing. [13, [1M| (Sec IVHl ) A nonzero result can be found in refs. l23l. [2^. With an alternative definition of spin 
currents, a nonzero result is also found in ref. fl^. 

For other odd N the SHE is nonzero [ISj, JJj] , assumes a universal value in a specific limit [l^ , but depends otherwise 
on the range of the impurity potential and is therefore in general not universal, though it is independent of the spin- 
orbit splitting 6, the overall strength of the impurity potential and of the impurity concentration [l^ . We reproduce 
these results in sec. lVIl Additionally, we give explicit results on the polarization and spin currents for the components 
of the spin parallel to the plane and calculate equilibrium spin currents. 

Wc will also show that the above results on the SHE are based on the response not being perturbative in the 
spin-orbit coupling. If perturbativeness is assumed (sec. IVHf|) . the spin-precession term — seemingly the prerequisite 
for a spin Hall current — must be left out to linear order in spin-orbit coupling, suggesting that to lowest order the 
intrinsic spin Hall effect is zero for arbitrary winding N (sec. IIX[) . Possible alternative contributions to the SHE from 
electric field-induced corrections to the collision integral arc discussed in sec. [X] A contribution that we believe has 
not been discussed before in the Boltzmann approach turns out to be the only candidate when the precession term is 
absent. However, this contribution turns out to be of second order in spin-orbit splitting. 

Leaving out the precession term leads to some formal difficulties in the case = ±1. The Boltzmann equation 
becomes unsolvable. However, this can be remedied by including a small spin relaxation term (sec. lIXp . This could 
suggest that for = ±1, non-magnetic impurities are not enough for a consistent steady-state solution in the spin 
polarization in the case that the response is perturbative in the spin-orbit coupling. A comparison of spin densities 
and spin currents in both phase space and real space is given by table [H In this context we also make a comment on 
pseudospin currents in bilayer graphene. 

The spin-orbit corrections to the collision integral make the analytical treatment considerably more complicated. 
For the case |A^| = 1 but not for the case \N\ I this complication is necessary even for a qualitatively correct 
understanding of the vanishing SHE. In appendix |B] we discuss the failure of the relaxation time approximation, that 
in many other problems is useful for a simple qualitative understanding and as a starting point for deriving real 
space equations for example describing thermoelectric effects. In particular, we contrast with the derivation of the 
Dyakonov-Pcrcl spin relaxation mechanism. 



II. THE MODEL — INTRINSIC VERSUS EXTRINSIC 

For a semiclassical Boltzmann description (see e.g. [13, [U, [13, [11] ) one needs the Wigner transformed one-particle 
Hamiltonian. For the spin-orbit (SO) coupled electrons we are going to study, it is 

H{x,p,t)=a{k)+(T-b{k) + e(j){x,t) (1) 

with e < and k{x, p, t) = p — eA{x, t). We want to describe a 2d system with x and p chosen to lie in the x, y-plane. 
Throughout the paper h = 1. In absence of spin-orbit coupling the dispersion is given by a oc , typically a = k'^ /2m. 
For the Rashba intercation b := \b\ = Xk and b := b/b = 9, with A parametrizing the strength of the coupling. 
However, we want to consider an arbitrary odd-integer winding number iV in 6 = b{k)b{0) (with bx + iby = 6*^"+'^^). 
The energy bands are = a -\- sb with s = ± giving the sign of the spin along the spin quantization axis 6, i.e. 
<j ■ 6|6s) = s|6s). 

The total Hamiltonian i?tot = H + Hi^p also includes an impurity potential Vimp{x) = U{x — a;„) of charged, 
non-magnetic impurities at positions a;„. Including the spin orbit coupling experienced at impurities one has 

Himp ~ Virnp H" Aoxt"" ' ^ ^ VVJmp ■ (2) 

H and -ffimp are treated very differently in the Boltzmann approach. H enters to linear order in the kinetic equation, 
whereas i?imp is in the Keldysh machinery turned into an impurity averaged self-energy to appear to quadratic order 
in the collision integral. 

A spin-orbit coupling enters both through the intrinsic {i.e. band related) term cr ■ b and in the extrinsic {i.e. 
impurity related) term Aoxt"" • k x VVimp, and the consequences of the two are usually studied separately in the 
literature. This paper deals only with the intrinsic contribution to the spin Hall effect. 
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III. SEMICLASSICAL DESCRIPTION OF A SPIN-ORBIT COUPLED SYSTEM 

In a Boltzmann description of an electron system with spin, the spatial degrees of freedom are treated semiclassically, 
whereas the treatment of the spin remains quantum mechanical. The state of the system is given by the 2 x 2-matrix 
valued distribution function f„„i{x,p,t), here with the spin index a ='\z^ iz- It is related to the equal time density 
matrix Pcrcr'{xi,X2)\t2=ti = {^\i{x2.,ti)'^ a{xi,ti)) by a Wigner transformation (see e.g. refs. [13, HU, [H, [H). In 
absence of scattering one can derive the Boltzmann equation for / by applying Heisenberg's equation of motion on 
p{xi,X2)j then identifying ^2 = ^i: Wigner transforming the result and gradient expanding it to first order. The 
approximation to stop at first order in gradient expansion is the semiclassical approximation, which relies on the 
external perturbations, such as electromagnetic potentials, changing negligibly on length and time scales of the de 
Broglie wavelength Ab and time tb = Xb/vf. 

From the matrix elements of the distribution function / one extracts the densities and current densities of charge 
and spin. The matrix elements are most conveniently expressed in the decomposition / = l/o + cr^/^ = /q + cr • / in 
Pauli matrices (with p, = x,y, z). (Throughout the paper we use the convention of summation over repeated indices.) 
Furthermore, we find it convenient to decompose the vector / = /gb + fgc + f^z in its components along the basis 

vectors b{9), z and c{9) = z x b{d), analogous to the cylindrical basis vectors k{d) := k/k, z and 6{9) z x k{9). 

The charge density en and current density ej in phase space are derived from en = Tr (fdH/dcj)) and ej = 
-Ti-{fdH/dA), which yields 

n{x,k,t) := Tr/ = 2/o = n+ + n" 

Mx,k,t) := Tr(v,;/) = 2fodk,a + 2f ■ dk^b = n+v+ + n'y- + ^ ' 

with i = x,y. Here we introduced the velocity matrices :~ di^-H — d^^a + cr • 9fc.6. The spin-independent part 
of the velocity is dka =: Vq. The band velocities are := dk('^ = {bs\v\bs) = v'^k. The intra-band elements 
:= (6± |/|6±) = /o±/g give the density of each spin band s = ±. The intcr-band elements (b± |/|b=p) = ±i/c 
are important for the coherent treatment of spin and are, for example, present in the last term of ^ [4y|, containing 
the Zitterbewegung of the spin-orbit coupled electrons. 

The real space densities are obtained by integrating the phase space densities over momentum, e.g. 

3{x,t)=: J -f^j{x,k,t). (4) 

When not otherwise stated, densities are in this paper always assumed to be phase space densities. 

The spin density, i.e. the polarization, is given by s'' = f Tr (cr^/) = (with ^ = 1). There is not a unique way 
to define the spin current because spin polarization is not a conserved quantity. (For a proposal on a conserved spin 
current, see ref. [s^. For its implications on the SHE, see ref. [l^.) When band velocities coincide, i.e. = Vq, then 
it is clearly = f^v^. For the general case we choose the common definition 

if = jTra^{v,, /} = f^dk^a + fodk\ (5) 
(with {A, B} = AB + BA). The spin Hall effect is a real space current of z-component spins 

jtTyI /^■"o = -c^SH^ X E (6) 
(27rj^ e 

perpendicular to the an applied electric field E along the plane, ctsh is the spin Hall conductivity. 
The Boltzmann equation in matrix form is given by 

i[H, /] + drf + ^{v„ d^J} + eE,dkJ ~ e..,ei?,^{v„ 9^^./} = J[f] (7) 

where the matrix-valued functional is the collision integral. In components it reads (from now on the charge e in 
eE and eB is absorbed into the fields) 

dtfo + dxJodk,a + dxj ■ dk,b + Eidkjo + ezijBzidkJodk^a + dkj ■ dk^b) = Jo 
2fxb + dtf + dxJdkM + dxjodk^b + E^J + ez,,Bz{dkJdk,a + dkjodk.b) = J (8) 

However, by virtue of definition ([5]) the equations ([8]) can be compactly written as 

dtn + d^-j + dk- [nE + j x B) = 2 Jo 
2(8 X br + dtsf' + ■ + dk ■ {s^E + j" xB) ^ J^. (9) 
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The left-hand side of the first equation is the same as for charged, spinless particles in an electromagnetic field. 
Apart from the spin-precession term, the second equation is of similar form. This is what one would expect since the 
electromagnetic field does not interact with the spin in the considered model but only with the charge that the spin 
sits on. The spin enters in a non-trivial way only through the precession term and through the collision integral. 



IV. DERIVATION OF THE COLLISION INTEGRAL IN THE KELDYSH FORMALISM 



The presence of, for example, two-body interactions or disorder averaged impurity interaction is in the Boltzmann 
approach described by the collision integral J . It is assumed that one is in the kinetic regime, where the de Brogiie 
wavelength A-r = l/fep is much shorter than the scattering length £. Like ref. we use the Kcldysh formalism (see 
e.g. refs. IsoL | 3ll. [32. fssi). but among other general methods we can mention the Nonequilibrium statistical operator 
formalism [30l|. For disorder averaged impurities, see also the compact derivation in ref. [2l|. 

The Keldysh derivation of the semiclassical equations starts with relating / to the Wigner transformed Keldysh 
Green's function 

f{x,p,t) = ^+ J ^^G'^{x,p,t,n). (10) 

The equation of motion for is given by the Dyson equation. Integrating the equation over the frequency results 
in the semiclassical Boltzmann equation with the collision term. This should be contrasted with the quasiclas steal 
Boltzmann approach (see e.g. ref. [ssh where the integration is instead performed over |fc| to obtain a distribution 
function f{x,p, t, fl) . This is the approach for example in refs. [H and [lo. 

The Keldysh equation derived from the Dyson equation can be written in two equivalent forms 

i = {idt-H-E)*G or i = G*{idt- H -Y.) . (11) 

The product involved here is the convolution product, the identity stands for 1 := IS {xi — X2) 5 (ti — 12) and quan- 
tities are written in Keldysh matrix space with 

G = A and idt-H -T.^\ w vA ■ (12) 



\^ G^J .^t ^ y Q idt~H -T. 

Each clement in these matrices is an infinite-dimensional matrix in real space indices, and in our case also a 2 x 2 
matrix in spin indices. 

The two equations in (jlip contain the same information. To derive kinetic equations one takes the difference of 
them, which for the Keldysh component yields 

{idt -H)*G^-G^* {idt -H) = * G^ - G^ * - G^ * + * G^ . (13) 

The right hand side is going to give the collision integral, which we only treat to lowest order in impurity concentration 
and impurity strength (first Born approximation), 

Sp = ni„p| ^|C/(|p-p'|)|2gO,, (14) 

with G^ = G(S = 0). (E is diagonal in momentum space due to the disorder averaging of the impurity interaction.) 
According to self-consistent Born-approximation we replace G^^ by G^. A crucial approximation comes with choosing 
the generalized Kadanoff-Baym Ansatz [sHi 

G"^ ^i{G^ *h~h*G^) + ... (15) 

where h := G^\t2=ti, which generalizes the quasiparticlc approximation G^ = h{x,p,t)S{il — e) for spinless electrons. 
The approximation can be considered to be an expansion in relaxation times of the system and corresponds to a 
Markov approximation. 

The convolution product takes after Wigner transformation A{xi, ii , a;2, ^2) ^ A{x,p, t, fl) the form A*B ~ Ae^^ B 
with the Poisson-bracket-like gradient T>. In a gauge invariant treatment valid when the electromagnetic fields are 
weak and vary slowly (see e.g. p. 344, Vol. 1 in ref. Isol ). one introduces k{p,x,t) = p — A and uj{Q.,x,t) = fl — (p 
and lets {x,k,t,uj} become the new set of independent variables {i.e. dx^k = 0). This changes the gradient into 

2? = dxidki- dkidxi+ d^dt- d t d ^ + Ei{d ^ d - d k, d ^) + e^JlBi^kj d k^ (16) 
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with X dY := {dX)Y and XdY := X{dY). Gradient expanding the left hand side of ([1^ to first order and 
integrating over the frequency yields the left-hand side of ([7]). The right hand side, the collision part, is usually taken 
to zeroth order in gradient expansion. Inserting (fT5|) into ([T^ . and using that combinations such as / duG^[. . .]G^ 
vanish, one arrives at the coUsion integral (with VF^^' :— 27rnii„p|C/(|fc — '^'DP) 

= -L^I ^(GfA/G".^ + G^? A/G°/) + . . . (17) 

where in the last row only terms of second order in the interaction strength were kept. The shorthand notations 
A/ = f{k,x,t) — f{k',x,t) and / ^^^2 J^' were introduced. 

In the expression (|17p we also need the retarded and advanced components. For this one should take the sum of 
the two equations pT|) . which for G^(E = 0) after Wigner transformation leads to 

2 = {u;+- i?°)e*^G™ + G°%i'^{Lu+ - i/O) , (18) 

with uj~^ := u! + irj (to take care of the boundary conditions provided by the imaginary part of when E ^ 0) and 
with iJ" := H — (j). To zeroth order in gradient expansion it is solved by 



Si,. „ 1 



G""=>^7T+^ 5g„:=;^(l + <T..b) (19) 



where S is the spin projection operator. With this zeroth order G°^, the last line in (IT71) can be reformulated to be 
identical to the result derived in ref. |2l| and essentially also to the one derived in ref. llOl . 

To the knowledge of the author, one finds in the literature always the zeroth order solution for G°^, maybe because 
in the common cases like spinless electrons or Zceman-coupled electrons (with a constant magnetic field) the first 
order contribution to (and thus to G°^) vanishes. However, for a spin-orbit coupled system it does not. We find 
the first order contribution 

„„nT} Eab — Bz(uj — a)dkb 

= -^^SE(-^-^-^^)^^ (20) 

s 

with := E-0 in the polar decomposition E = E^k + EgO. In section|X]wc show how the contribution modifies 
the collision integral. We will also investigate, whether the corresponding correction could be an alternative source 
to the SHE not relying on the precession term.[47| 



V. COLLISION INTEGRAL TO LINEAR ORDER IN SPIN-ORBIT COUPLING 



In this section the Boltzmann equation is expanded to first order in the spin-orbit coupling, as done in refs. Il9l andl2ll. 
hence assuming b{k-p) <C ep- The subscript F indicates the value of the corresponding quantity at the Fermi surface 
determined by a = ep, where ep = ^ at low temperatures fceT <C ep. 

The collision integral is taken to the habitual zeroth order in gradient expansion, meaning ()17p with (|19p . To 
slim down the often lengthy expression for collision integrals, some more shorthand notation is introduced, x' means 
that the quantity x depends on primed variables such as fe', s' etc, whereas x correspondingly depends on k, s. For 
example S" = ^(1 + cr ■ s'bfc')- Also, Ax x — x', for example Ae = — e^, or A{sb) ^ sb — s'b' . 



Inserting into ([T7]) gives 



Jo = - I W,,k' ^E^(Ae) 

ss' 



1 + ss'b b' ^ sb + s'b' ^ ^ 
A/o + A/ 



i±i£^A/ + f^±i^A/o 



^0 



X . 



(21) 
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with the matrix B;,^.' := b{b')"^ + b'{b)^ — b ■ b'.[4^ Not written out in equation (|2T|) are the principal part terms 



X 



ss'b X b' , , A(sb) 

(A/ - <tA/o) - a ■ X A/ 



(22) 



which are not considered to be a part of the elastic collision integral but to be related to renormalization corrections. 
Such terms can usually left out if one is interested in the interaction only to lowest order (however, see ref. IstI ) and 
it is beyond the scope of the present paper to discuss them. See ref. [s^ for a lucid treatment on how to interpret and 
handle them in the case of spinless electrons interacting through a two-body interaction. See also the derivation of 
Bloch spin relaxation equations in ref. [sol . 

The delta functions (5(e| — e^/) connecting Fermi surfaces at different |fc| make it difficult to find an analytical 
solution. However, in the considered limit 5(fcF) <C ep one can use the expansion 



(5(Ae) = J(Aa) + A(s6) (5'(Aa) + 0{\^) 



(23) 



after which one is left with the spin-independent delta function 5{ak 
the collision integral reads 



a/j/) implying — \k\. With this expansion 



Jo = - f Wkk' [5{Aa)Afo+'6'iAa)Ab-Af) 

Jk' 

J = - I W^k' [S{Aa)Af + S'iAa)AbAfo] . 



(24) 

Jk' 

The terms with S'{Aa) = --da'S{Aa) are made sense of by integration by parts. 

If / = 0{X) {i.e. if the polarization vanishes as A — > 0), the term indicated as of order 0{X^) can be neglected to 
linear order in A. This resulting collision integral is essentially the one found in ref. [2l] and is similar to the one in 
ref. [iqI. [iot For a physical interpretation of spin-orbit coupling dependent contributions, see ref. [l9l. 

In section|X]we discuss corrections to the collision integral when one goes beyond zeroth order in gradient expansion. 



VI. SOLVING THE BOLTZMANN EQUATION 



We now set out to solve the uniform, steady state Boltzmann equation 

2a- f xb + E-dkf = - I VKfcfc' [5{Aa)Af + 5'{Aa)Ab-(TAh] (25) 

Jk' 

without implementing the assumption of perturbativeness in A, to be discussed in sec I Villi The distribution function 
/ = + f^^^ is linearized only in the electric field and not in the spin-orbit coupling. 

The distribution = /o ± of respective spin band is in equilibrium given by the Fermi-Dirac (FD) distribution. 
The equilibrium distribution is therefore given by 

fl'' ± /^^ = /fd (6± - ^l) i.e. r - E ^6±/fd (^^ - m) , (26) 

± 

with vanishing inter-band elements = f^"^ = 0. Time reversal symmetry requires the real space equilibrium 
polarization = f^^{k) to be zero, which also follows trivially from the vanishing angular part of the integral. 
The real space spin current, on the other hand, need not vanish since it is even under time reversal symmetry. For 
A^l 7^ 1 the real space spin current is trivially zero, but for N — ±1 we find jy = =Fj^ = —w?\^/2'n + O(A^) for a 
quadratic dispersion a = /2ra. 

The equilibrium distribution is from now on taken to linear order in A, i.e. 

C = E± ^^^^ = /fd (a - M) + 0{X') ^ f,'^ = /fd ,2^. 

^cq ^ ^ ^^^^^^ _ ^ n = bdahu 

where /fd = /fd (^ ^ m) fi'om now on. From (j27|) one sees that a small spin-orbit coupling does not change the charge 
density but induces a small polarization at the Fermi surface. The distribution f^"^ = /fd + f • bda fpD satisfies 
for E^O. 
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The charge part of ([25)) does not depend on the polarization. When a uniform, static electric field is applied, one 

finds the usual solution /q^^ = —TtiE^dk/pD where Ttr is the transport relaxation time. Gathering the known terms 
on the left-hand side one can write the polarization part of the equation as 



E ■ dkf 



oq 



k' 



k' 



It turns out (see the appendix) that the equation can be written in the form 



F{E.b - Esc) + G{Erb + Esc) = 26 x / 



{E) 



KAf 



{E) 



\k'—k 1 



(28) 



(29) 



with the shorthand Jg, := / The functions F and G depend on k and are proportional to A. They do not depend 
on 9, on E, on the impurity concentration njmp or on the overall strength of the impurity potential U. However, they 
depend on the range of the potential through dimensionless fractions of the Fourier components of the functions 



K{k,A9) 
k{k,A9) 



in 

D{a)kvo[da'Wkk']k'=k ^ 



imA.9 £> 



(30) 



[5(j | Particularly, for = 1 we find for a general dispersion a cx k'^ that 

9a /fD 



F ^ \ 

G = 



((C-2)r,7i+Cri 



-1 

2 



2TtrT]^2 ^01 



'01 



'12 



Ttr + 



kdkdaf] 



FD 



'l+TtrTia') 



(31) 



with = Kq — Ki and introducing shorthands := A'l — K2, f^^ := Kq — Ki and ~ ^ 1^2- We find that 
G = also when 6 is not proportional to k. The vanishing of G is going to imply the vanishing of the spin Hall effect 
for the Rashba coupling. 

The combination EjJ) =F E^c has a winding ± 1. Particularly, for the Rashba case b ^ {i.e. N ^ 1) one has 



Ej^b-EgC= [b{k)' - c(6/) 



-sin2A6l cos2A6' 
cos 2A6' sin 2A6' 



Ej^b + E^c^ 



-1 

1 



(32) 



i.e. the left-hand side of ([25)1 has an angularly independent term. 

The solution can be found by Fourier decomposition in the basis {6, c, 2} 




(33) 



where the Fourier coefficients {/;,„, fen, fzn} of course only depend on k and not on 9. With £ := E^ + iEy one has 
2{Ej^b ± Egc) = e'®£*(6 ± ic) + e^*^£(6 =F ic) and therefore the left-hand side of ([29)1 can be written as 



-€'"£* {bcz)\tG-tF 



F + G 




c.c. , 



(34) 



which contains only the n = 1 Fourier component and the complex conjugate n = —1 component. This is going 
to imply that /g^^, fen, fzn ~ for \n\ ^ 1. (Choosing a cartesian basis in ([33]), in contrast, couples the equation 
for component n with the components n ± N.) This fact has some direct implications for the electric field induced 

contributions to the real space densities. For example, the real space density of z-spins fz is trivially zero. For 

( E) ( E) 

\N\ ^ 1 the in-plane components fx and fy also vanish trivially in real space, whereas for \N\ = 1 they can be 
nonzero. The real space spin Hall current j^' can be nonzero for all N, whereas the contribution to and vanishes 
trivially for all N. 

On the right-hand side of ([25)) we find 



6 (4,^ Jg, K{1 - COS NA9 cos nA9) + if an Jg, K sin NA9 sin nA9 
c [f&n Jg,K{l- COS NA9 cos nA9) - if-^^^ Jg, K sin NA9 sin nA9 
z fzn Jg, K{1 - cos nA9) 
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Here, b' = b cos N A9 — csm N Ad was used and terms odd in AO were left out, using that K is even in A6. Including 
the precession term 26(i/c — c/^) the equation for the 72=1 Fourier components of /(^' becomes 



£* 

T 





+ G \ 




iG 


-iF 






/ 






(35) 



where 



:= / -cos7VA6lcosA6i) = A'o - (i^w-i + i^w+i)/2 
Je' 

tT^^I := I KsmNA9sinA9 = {Kn-1- Kn+i)/2 (36) 



Je' 

(possibly negative) were introduced. 

With the A-dcpendent b, the matrix and hence the solution will be inhomogenous in A. For jA^j 7^ 1 the solution 
nonetheless goes to zero when A — > 0. For point-like impurities {t^J^ = '''tr^ = ^'^"0 a-nd rTjj^ = 0) we find 



4fc2A2 + if2 7^-^"^^ 27r(4fc|A2 + if2)V ^v^- 

for the real space spin current. As in ref. [l^ one recovers in the clean limit r^^^ = Kq <^ fcpA a universal spin Hall 
conductivity cr^^ = whereas ctsh decreases to zero in the opposite {dirty) limit. 

For the case N = ±1 one has for arbitrary impurity range that Tcos = i^sin 2ro2- The determinant Ab'^r^^ + 
'''tr ^ ('''cos ^ '''sTn) "f the matrix in ((35|) becomes singular at 6 = 0. For iV = 1 we find the solution 

h^ = -^*(-02(F + G) + ^ 



hi = 



G 



4Tt,62 



/.I = -^^*Yb- ^^^^ 

The z-component of the real space spin density must vanish trivially as noted after equation p4p . The in-plane 
spin components, resulting in a electric-field- induced polarization in real space 

s = ^(b/g + c/e) = -zxE j daD (^to2{F + G) + . (39) 

Remember that the contribution from the equilibrium polarization is zero in real space. 

According to (|3ip . G = for N = 1 for arbitrary dispersion and arbitrary range of impurities according. Thus, 
fz~0 for a Rashba coupling, implying a zero spin Hall current. A similar analysis applies to = —1, e.g. for the 
linear Drcsselhaus coupling, since a model of winding ~N is related to one of winding N through reflection (say along 
the X-axis). 

If for the Rashba coupling there was a nonzero spin Hall effect, i.e. G ^ 0, then jz and consequently the spin Hall 
current would be indpendent of A (since in general, G oc A). Furthermore, the in-plane components and /j, of the 
polarization would diverge as A~^ as A ^ 0, i.e. one would not recover an unpolarized distribution when sending the 
spin-orbit coupling to zero. However, in the derivation of the Boltzmann equation there were only assumptions of A 
and / being small enough, no assumptions of them not being too small. Therefore, to the extent that a diverging 
polarization in such a case is an unthinkable result, the vanishing of the SHE for the Rashba case is a natural 
implication. 



VII. GENERAL ARGUMENTS FOR A VANISHING SHE 



The vanishing of the SHE for the Rashba case has been found by numerous previous studies (see the introduction). 
To the knowledge of the author it has not been related to the finiteness of the in-plane polarization like done in the 
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argument above. In linear response studies general arguments have been given in the case of quadratic dispersion 
a (X [13, nil, where it is noted from [ay,X{p x tr)^] = iiXpyOz that for Heisenberg operators 

^ ~i[&y,H + Vin,p] 2Xpya, ^ 2\mfy (40) 

for non-magnetic impurities. The steady state condition {■^o'y) = forces the spin Hall current (j.^) to be zero. A 
similar argument applies to . 

The analog to this argument in the Boltzmann approach can be established for the real space densities. In phase 
space, (s X b)^ = bxS^ = jybx/vQy. Only for = ±1 and for b/vo independent of k is it possible for bx/voy to 
be a constant, here ^Am. According to ^ one has dts^ = — 2(s x b)^ — E ■ dkS^ + in the uniform case. In 
the integration over momentum the last two terms vanish, resulting in dt-s^ ~ ±2Xmjy in real space for A^ = ±1. 
Likewise, dts^ = —2Xmj^. In a steady state, dts = implies = 0. 

In phase space, on the other hand, the steady state condition = dfsy = — 2(s x b)^ — E ■ dkS^ ■ ■ ■ does not imply a 
vanishing jy. We could have had G as long as / da DvoG/b = guaranteed the vanishing in real space. However, 
in the previous section we find that = also in phase space. 



VIII. ASSUMPTION OF LINEAR RESPONSE IN THE SPIN-ORBIT COUPLING 



We are now going to study the implications of a basic assumption, namely that the steady-state response of a spin- 
orbit coupled system in an electric field is perturbative in the small parameter A {vp is kept constant), analogously 
to the usual assumption of linear response in the electric field strength E. For the latter assumption it means that 
one can expand the solution f^'^ + f^^^ + f^^ in nonnegative powers in the electric field and solve the equation 

iterativcly by solving equations that are homogenous in orders of E, see e.g. equation (j28p . With the assumption 
now to be studied, the same is expected to apply also for the spin-orbit coupling parameter A {v-p is kept constant), 
wherefore the solution can be expanded in powers of E and A 

/= E /^'"''"^ (41) 

m,n>0 

assuming that the solution is analytic both in A and E. 

The equations are going to be solved order by order in both parameters as illustrated in (|45p . In the expansion of 
the Boltzmann equation in a small spin-orbit coupling, there is in this case no ambiguity about which terms to include 
in a given equation. This should be contrasted with the previous section, where we chose to include the precession 
term 2b x f — a term of order O(A^) — though we discarded other terms of the same order, for example in ([24)) . The 
reproduced results seem to be based on hand-picking terms with physical insight. 



IX. BOLTZMANN EQUATION WITHOUT THE PRECESSION TERM 

Due to the assumption in sec. IVIIIl not only the equation but also the solution is linearized in the spin-orbit coupling. 
Together with the habitual linearization in the electric field this means that we consider 

./ = ./^°^ + ./^^' + Z^'"' + .f^^-""^ (42) 

where the superscripts denote the order in A and E, respectively {e.g. Z'^' oc X^E^). The equilibrium distribution is 
n = /^°^ + /^^\ where /(o) = I/fd- 

The left-hand side of the Boltzmann equation ([5]) is in the static, uniform case 

%=E-dk.f + 2a-TTb (43) 
at 

where the precession term must be neglected to order A. However, leaving out the precession term leads to some 
formal trouble in the case A^ = ±1. The derivative E ■ dk = e'^^{Ex — iEy){dk + ik"^de) + c.c. comes with a winding 
number ±1. For a spin-orbit coupled system the equilibrium polarization /°'^ = f{E = 0) has a also winding number, 
here A^. Thus, the combination E ■ dkf comes with terms of winding A^ ± 1. For A^ = 1 {e.g. Rashba) or A^ = — 1 
{e.g. linear Dresselhaus) the left hand side therefore contains terms without angular dependence, as seen in ((32|) . 
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Such terms cannot be matched by the colhsion integral, essentially because an equation like 1 = / d6' {f{6) — f{0')) 
has no solution. (A collision integral cannot be a source/drain.) 

This lack of solution might be related to the discarding of the principal part terms X in ([2T|l . but we have not been 
able to investigate this. As a simple remedy, we introduce instead a small spin relaxation term, which could come 
from spin relaxation processes not related to the spin-orbit coupling. The Boltzmann equation then reads 



rr-Ts'f + E-dkf = [ I^fcfe' [S{Aa)Af + S'{Aa)a-AbAfo] 

Jk' 



(44) 



where it is assumed that t^^ is much smaller than t^~^, is independent of A and is a number though in general it 
could be a matrix. 

We now solve the Boltzmann equation (j44[) order by order in A and E as shown in equation (|45p 



X°E^ 
Ai^i 



= -/fe'M^fcfc'^(Aa)A/ 



(0) 



(T ■ n 



'V*^^ = -L'Wkk' SiAa)AfW+S'{Aa)a-AbAh 



(0) 



E ■ = - h, W^^, [<5(Aa)A/(^^^) + S'{Aa)a- ■ A6A/f ' 



(45) 



Note that terms known to be zero were left out.[5l| We are out to find /''^^^■', which is of order X^E^ and therefore 
the first contribution that incorporates the combined effects of an electric field and spin-orbit coupling, [s^l However, 
to solve the X^E^ equation one needs to solve the previous equations to find out /^^^ and f^^\ 

The X^E*^ equation is consistent with /'^^^ = 1/fd from ((27|) . The spin-relaxation term in the X^E^ equation 
decreases the equilibrium polarization in (|?7|) into ~ 70- ■ bdafpi) with 7 := (1 -I- r^^rtr)"^ ~ 1. The X^E^ 



equation is solved by /'^•' = — l7tr-£'£.(^fc/FD- So finally at the X^E^ equation one knows already the terms E ■ d^f^^^ 



and /^/ W^,^'6'{Aa)Ab ■ (tAJq ' . Collecting these contributions on the left-hand side we get an equation of the same 
form as (j29|) , but with the spin precession term replaced by the spin relaxation term, and with F and G modified due 
to the factor 7 in £^ • dkf^'^^ ■ The equation for the n = I Fourier coefficients turns into 





+ G \ 




iG 


-iF 






/ 





cos 

.-1 













(hi 




fcl 




{hi 



(46) 



For = ±1 one has r. 



unless Tg 7^ 0. 



cos 

(With 



for a general impurity potential, leading the determinant of the matrix to be zero 
there is either no solution when G 7^ or multiple solutions when G = 0.) For |A^| 7^ 1 
one has Tcos ^ Tgin for realisitc impurity potentials, wherefore t^^ 7^ is not needed. 

Important to retain is that the matrix does not depend on A. The polarization is therefore proportional to 

A. Clear is also that /z = and the spin Hall current j'^ are zero — for an arbitrary spin-orbit coupling. 

For A^ = ±1 we can adopt the derivation in sec. IVIII to show that dtSy = ±mXjy — Tg^Sy in real space. For 
our steady state case we therefore find cx = in real space for A^ = ±1. Likewise, = 0. For |A^| 7^ 1 
the real space polarization is trivially zero (see under equation p4p '). Summarizing we have that for no A^ in the 
perturbative case does the electric field alter any real space densities to lowest order in A. The only nonzero real space 
densities are the equilibrium spin currents and for |A^|=:1. From f^^ = /fd + 7c • bdafpi) for A^ = ±1 wc find 
= TJ^ = 4^(7 - 1)^fA:f. 

Table H] gives a summary and offers a comparison of the different cases consider in this paper. In this paper we do 
not consider the case A^ even, but we note that for |A^| = 2 an electric- field-induced contribution to the real space 
spin currents and is not forced to be zero by simple angular considerations. (This is analogous to s^*-^^ and 
being allowed to be nonzero only for |A| = 1.) The case A^ = 2 is relevant for pseudospin currents for one valley 
in bilayer graphene. Note that this is not the SHE pseudospin current discussed in ref. H^. Actually, with the 
definition (O a nonzero (pseudo)spin current is not possible when a = 0, the latter being the case in single and 
double layer graphene. 



X. ELECTRIC FIELD INDUCED CORRECTIONS TO THE COLLISION INTEGRAL 



The precession term seemed so far like the only term that could involve the fz component and lead to nonzero spin 
current . In this section we arc going to see that the electric field E modifies the collision integral in a way that 
involves the fz component. However, the corrections turns out to be of order A^. 
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Nonperturbative 





\N\ = 1 


fiE) 


NZ/ 0/0 NZ/ NZ^' 0/ 
NZ/ NZ^) O^V NZ/ O^V O"*' 


\N\ / 1 


ycq 
fiE) 


NZ/ 0/0 NZ/ 0/0 
NZ/ NZ/ NZ/ 0^> NZ/ NZ 
Perturb at ive 


|Af| = 1 


jcq 
fiE) 


NZ/ 0/0 NZ/ NZ^) 0/ 
NZ/ 0/ NZ/ 0/0 


|Af| ^ 1 


joq 
fiE) 


NZ/ 0/0 NZ/ 0/0 
NZ/ 0/0 NZ/O'^) 0/0 



TABLE I: Summary of spin densities and spin currents in the nonperturbative and the perturbative cases for the cases = ±1 
and for other odd N. The entries present the (phase space)/(real space) quantities, respectively. NZ stands for a (typically) 
nonzero result whereas uncommented zeros stand for trivially vanishing components {e.g. due to angular considerations). 
Nontrivial cases are commented. 1) j^,jy = C'(A^), see under (|26|) . 2) See p9[) . 3) Since G = 0. 4) General arguments in 
sec. IVIII 5) Vanishing for N odd. Could, however, be nonzero for |A| = 2. 6) See end of sec. IIXI Relies on the introduced spin 
relaxation term. Here j%,jy — 0{\) 

In trying to incorporate in our Boltzmann equation all terms present to first order in our parameters we so far 
left out terms by gradient expanding the right-hand side of the Dyson equation (fT3|) for only to zeroth order 
and not to first order in the electric field. First order corrections have been accounted for example in the case of 
electron-phonon renormalization of the ac conductivity [4l| (see also [s^l) and have also been discussed in the SHE 
and AHE literature (see e.g. refs. [H orl43). However, a derivative such as 2E ■ dkSf^^ = sk~^Egdgb ■ cr, which occurs 
in the derivation of these corrections, gives a vector that remains in the plane. The collision integral becomes more 
complicated but does not involve the fz component. 

The contribution discussed in this paper comes about in a slightly subtler way, and has to the knowledge of the 
author not been discussed previously. It comes from gradient expanding the equation of motion also for the retarded 
Green's function G^^ , which results in the correction 5G'^^ given by (pOj) . With = one obtains 

SJ ^ -J^^^J ^(JGf A/GO^ + Gf A/<5G°^ ^ G^?A/<5Gf + SG'^AfGl^) = 

( a,AS' + S'Afa,)Eg{da- sb-^)SiAe)+ 
+ (a, AS + SAfa,)Eg,ida' - s'b'-')6iAe) _ 

= <T • 1 5"iAa) [W^^- [zM ■ Af + MA/,]) (47) 

where for the last line the expansion (j23p was used and where 

M(fc, k') := k-^Egb' + k''^Eg,b . (48) 

Note that 5J^q = 0, i.e. there is no contribution to the charge part of the equation, only to the polarization part. 
Note particularly that 5Jz 7^ Oi which would give a nontrivial equation for However, since M oc \E and / cx A, 
this correction contributes to order \^E as announced, wherefore the result fz~0 in last section is not changed to 
lowest order. 

One can ask if this correction could change the result ~ for the Rashba case in the non-perturbative scenario 
of sec. I VII where homogenity in A was not an issue. The correction would then enter as an imaginary element of order 
X^E replacing the zero in the bottom of the vector with F and G in ([55]) . However, since t~^1 — t^^ = for A^ = 1, 
the inverse matrix has a zero zz-element. Therefore, the correction cannot contribute to fzi in the Rashba case. In 
this aspect the contributions discussed in the beginning of this section, on the other hand, could contribute, but we 
have not investigated them in a systematic way. The real space spin current would in any case remain zero, at 
least for a quadratic dispersion, due to the arguments in sec. IVIII 
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XI. CONCLUSIONS 



This paper studied the intrinsic contribution to the spin Hall current in a spin-orbit coupled 2DEG by deriving 
a Boltzmann equation in the Keldysh formalism and solving it in the uniform steady state case. The vector b 
determining the spin-orbit coupling was assumed to be of the form b = b{k)b{9) with b = Xk and bx + iby oc ^^^o+^Ne ^ 
We reproduced the common result that spin Hall effect vanishes for N ~ ±1 {e.g. for a Rashba coupling) but not for 
other N . We were able to give a new perspective on this vanishing by pointing out that a nonzero result leads the 
in-plane components of the polarization to diverge when A ^ 0. 

The mentioned treatment does not assume the response to be perturbative in A. We therefore found it interesting 
to study the implications of assuming the response to be perturbative not only in the electric field but also in A. 
The precession term — previously the prerequisite for the spin Hall current — must then be left out to first order in 
spin-orbit splitting. The out-of-plane polarization becomes trivially zero and there seems to be zero SHE for any 
winding. We saw also that all other real space densities have zero electric field induced contributions. 

Leaving out the spin-prcccssion term gives a Boltzmann equation for in-plane polarization {fx, fy) that is unsolvable 
for N = ±1. The unsolvability might be related to the left out principal parts, the inclusion of which would have 
been beyond the scope of the present study. As an ad hoc remedy the precession term got replaced by a small spin 
relaxation term. This could suggest that if the response is perturbative in A, then non-magnetic impurities are not 
enough for the existence of a steady state solution for the polarization. 

To cover all contributions to first order in electric field, in spin-orbit splitting and in impurity strength and con- 
centration, we considered corrections to the collision integral that come from going to first order in electric field in 
the gradient expansion of the self-energy side of the Dyson equation. One of the corrections, to our knowledge not 
discussed before, actually involves the fz component. However, this contribution is of order A^. Thus, the vanishing 
of the spin Hall current to lowest order in A in the perturbative case is not changed by these corrections. 

Finally, the paper includes a detailed discussion of why a relaxation time approximation fails and a comment on 
pseudospin currents in bilayer graphene. 

Aknowledgements. The author wishes to aknowledge discussions with M. Liiffe, D. Culcer, F. Gethmann, 
A. G. Mal'shukov, K. Morawetz, T. Nunner, P. Schwab, G. Vignale, F. von Oppcn and R. Winkler. This work 
was supported by the Swedish Research Council. The author also wishes to acknowledge a visitor grant of the Max 
Planck Institute for the Physics of Complex Systems. 

APPENDIX A: LEFT HAND SIDE OF EQ. (f29]) 



For b — Xk and b^ + iby = e'^o+'^^ and / — ^bda /fd one has 

1 

k 



E-dkf^ {E-^dk + lEgdg)jXkbda fFB - jX{Ej^b + NEgc) dafFB + iXk Ej^bdkdafFD ■ (Al) 



The term _|^/ W/^i^'S' (Aa)Ab ■ ctA/q^^ is here for brevity only evaluated for a point-like impurity potential, i.e. 
Wj^y = W constant, and constant density of states D{a) = TO/27r. Hence 

r^;! = / 5{Aa)W^^,{l - cos(A0)) = DW (A2) 

Jk' 

With 5'(Aa) = da5{Aa) = ~dai5{Aa) a partial integration gives 

j WS'iAaJAbAf^'^^ ^ j da''^5{Aa)da'{DWAbAff'^)^ 

da'SiAa) kdadkfFD / —E^,Ab + da'{k')dk.fFD / i^^'AE^, = 

= -S\N\,l^iEi,b + NEgC)ikdadkfFD + dafFD) = 

= -d\N\.l^{Ef^b+NEgc){kdkdafFD+CdafFD) (A3) 

where it was used that b' = bcosNAO — csiiiNAd and k' ~ k cos Ad — 0sinA0. In the last line it was used that 
dadk ~ dai^da) = {daVo)da + ^od^ = ^^^a + i9fe9a for vq OC k^~^ . For a non-constant D {i.e. for C 7^ 2) the result 
in (jA3p is modified. Note also that it is only for point-like impurities that the contribution (|A3[) vanishes for |A^| ^ 1. 
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Adding up and (|X3|) one obtains for TV = 1 that F = ^X{C - '2)dafFr> + ^^ikdudahi) and G = \{i - 

l)(9a/FD + ■^kdkdafFi})- For jA^I 7^ 1 the contribution (jA3p vanishes and one gets for example G — F ^ NX'^dafp'o, 
needed for the result ([57)1 . 

APPENDIX B: FAILURE OF A RELAXATION TIME APPROXIMATION 

A useful approximation for the collision integral found in standard applications of the Boltzmann equation is the 
relaxation time approximation (RTA) 

Jlh] = - I 5{^e)Wk^>{h - h') -— , (Bl) 

Jk' T 

where the relaxation time r depends only on the absolute value k of the momentum. It expresses that the role of the 
collision integral is to relax the deviation from equilibrium Sfk '■= fk — fk^- For M^^fc' angularly independent (i.e. 
momentum independent) it can be derived from the collision integral to the left in (|Bip . Thus RTA should be a good 
approximation for weakly momentum dependent impurity potentials. 

A prerequisite for the RTA is that it is consistent with the conserved quantities. In particular, the momentum 
integral of the Boltzmann equation should give the continuity equation dtn + dx ■ j = expressing the conservation of 
particle number. The right hand side (given W^^' = W^'f,) indeed vanishes identically J^^ J[fk] = 0, expressing that 
the collision term cannot in real space act as a source (or drain) of particles. For the RTA the particle conservation 
is not automatic, since Sfk/r would be nonzero if Sf contained an angularly independent component. However, in 
typical applications the deviation Sf = ^„ (5/„ exp(m0) contains to lowest order only angularly dependent terms. 

The RTA can be implemented also in the case of spin in the same way as in (jBip . or more generally by letting r^^ 
be a matrix acting on 8f to allow for different relaxation times for different spin components of /. The RTA has been 
useful for example for the derivation of the Dyakonov-Perel spin relaxation mechanism (DPSR) [i^ (see also refs. H 
and [45I). found by deriving for initially polarized electrons the Bloch equations from the uniform, time-dependent 
Boltzmann equation for a Rashba-type spin-orbit interaction. In that problem the electric field is absent. 

The general collision integral pT|) still satisfies = for all components, expressing that not only particle 

number but also spin is conserved in collisions with nonmagnetic impurities. Thus, a sensible RTA must still satisfy 
/fc ^fk/T = 0. That is to say 5f cannot contain an angularly constant term. In the derivation of the DPSR this is 
satisfied although n = components are present and essential. The RTA is only needed in the part of the Boltzmann 
equation that is of first power in A, and Z^^-* is by simple inspection seen to only contain the angular components 
n = ±1. 

We note that in the SHE setup studied in this paper the situation is very different compared to the DPSR problem. 
It turns out that a RTA always implies a nonzero SHE. For = ±1 we therefore do not even qualitatively reproduce 
the correct spin current. For |A^| ^ 1 the SHE is on the other hand nonzero and can be captured by a RTA. In 
the A = ±1 case the nonzero SHE is intimately related to the RTA failing to conserve spin. With the simple RTA 
J[f] —Sf/r the real space equations for the Rashba/Dresselhaus case in section IVHl get modified to 

dtsy ^±2mXj^~Ssy/T. (B2) 

A steady state does no longer imply jy = but instead jy = dzSs^ /2mXT. Thus, if a polarization in real space is 
induced by the electric field, then there is a nonzero SHE. The result in (|39p shows that such a polarization is actually 
induced. Thus, the collision integral transformed into RTA form would enter as a nonzero spin source, simultaneously 
allowing for a nonzero SHE. Therefore, in contrast to the DPSR problem with a Rashba coupling, the RTA fails when 
applied to the SHE problem with a Rashba coupling. 

For N = ±1 in the perturbative case, a RTA and equation (jB2p do not lead to the same contradiction since 
= = in real space (see table |l|. In the rest of the section we give some further details on the nonperturbative 
case. 

In the SHE problem we look at a steady state situation and the expansion is done around the equilibrium distribution 
instead of an initial polarization as in the DPSR problem. The Boltzmann equation is expanded in powers of E rather 
than in powers of A. The lowest order deviation /'■^■' (from the equilibrium polarization f'^'^) will now contain 
angularly independent terms, causing the non- vanishing of J^, Z^^'' in (|39p . Contrast with ordinary steady problems 

with an electric field, where constant terms are not present in Z*^-^^ ~ tE ■ vdef°'^ since f'^'^ is independent of 9 and 
E ■ V comes with angular components n = ±1. In the SHE problem on the other hand f^"^ oc b comes with angular 
components iN, and in the Rashba/Dresselhaus case in particular N = ±1. In the Rashba case Z'^-* can therefore 
contain angular components n = ±1 ± 1 = —2, 0, 2 in the cartesian basis. (In the rotating basis given by {b, c, z} used 
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in this paper it translates into the Fourier components n = ±1. ) Particularly, the presence of the n = term means 
that a RTA would be inconsistent with the impurity scattering having to conserve spin. 

The other side of the coin — as displayed in (jB2p — is that the RTA allows for fzi^O and consequently a nonzero 
SHE. The analog of equation (j35p is with a general matrix with constant elements is given by 



£* 
T 



V / 



/ 



V 



bb 



'be 



-2b 



^-1 







1 


hi 







(B3) 



Note that now the left hand side derives only from the term E ■ dkf^^, whereas the left hand side of ([HS)) had an 
additional contribution from the term J^/ W/^i,' 6' {Aa)AbAfQ^^ as seen in (|28p . This latter term cannot be captured by 
a RTA. However, this term was essential for giving a left hand side in (j35|) the crucial structure that all the components 
of the vector are proportional to each other, i.e. linearly dependent, as followed with G = 0. This structure makes 
it possible to find a solution f^i = for suitable matrix elements in r. (In (|35p we had Tcos = iTsin-) In the 

f^^ are linearly independent. 



vector of the left hand side of (jB3p . on the other hand, the components dkf^ and k ^ 



There is no natural choice of matrix elements of that with such a left hand side can result in the cancellations 
needed to make f^i vanish. (The elements of are assumed to be at the most weakly dependent on k and should 
certainly not contain factors of Fermi-Dirac distributions.) For the same reasons also a simple collision integral like 
J ~ Wkk'^{^o)Af fails to reproduce the vanishing SHE because it gives a Boltzmann equation with the same 
left hand side as in (|B3p . 

For the case |A^| ^ 1 the term J^,, W;,^/(5'(Aa)A6A/g^^ vanishes for W constant (corresponding to t constant), 
as seen in appendix [X] The collision integral enters the Boltzmann equation only with the simple contribution 
/fc' ^kk'^{^'^)^f 1 which for W constant can be put on the RTA form. Therefore, the RTA can qualitatively reproduce 
the correct result. Neither can there be any angularly constant terms in f'^^'^ for \N\ ^ 1. Here the RTA works well, 
in contrast to the lA^I = 1 case. 
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